Novel Solution of Mercury Perihelion Shift 
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We present a novel solution ol the Mercury perihelion advance shift in the new gravity model. It 

is found that the non-relativistic reduction of the Dirac equation with the gravitational potential 

. , rrr, ^ GMiu ^, . 

produces the new gravitational potential oi V(r) = 1 . i his potential can explain 

the Mercury perihelion advance shift without any free parameters. Also, it can give rise to the 
cij— shift of the GPS satellite where the advance shift amounts to ( — ) ~ 3.4 x 10"^*^ which should 

\ ^ / th 

be compared to the recent observed value of ( — ) — 4.5 x 10^^". 
PACS numbers; 95.30.Sf,95.30.-k,11.10.-z 



I. INTRODUCTION 

In a recent paper, a new gravitational model is pro- 
posed which can well describe all the known properties 
of the gravity physics [U, In addition, this model 
predicts a new interaction between photon and gravity 
at the fourth order Feynman diagram, and this photon 
gravity interaction should play some important role in 
astrophysics. 

From this gravity model, we obtain, for the first time, 
a proper Dirac equation with the gravitational potential, 
and it is written as 



-iV ■ a 



(1.1) 



where G and Af denote the gravitational constant and 
the mass of the gravity center, respectively. This equa- 
tion can be reduced to the non-relativistic equation in 
terms of Foldy-Wouthuysen transformation, and after we 
make the classical limit, the new gravitational potential 
becomes 



V{r) = - 



GMr 



G^M^m 1 
2c2 ^ 



(1.2) 



This additional potential is special for the scalar po- 
tential, in contrast to the Coulomb potential which has 
no such effect from the non-relativistic reduction of the 
Dirac equation. 

In this paper, we show that the new gravitational po- 
tential can explain the advance shift of the Mercury per- 
ihelion to a remarkably good accuracy. Also, it gives rise 
to the advance w— shift of the GPS satelhte (—),^ ~ 

V uj I th 

3.4 X 10"^" which should be compared to the observed 
value of (^) -4.5 x lO^i" 131, and one sees that 

\ UJ / exp ^ J ' 

the theoretical evaluation achieves a remarkable agree- 
ment with experiment. Surprisingly, all the calculations 



are carried out without any free parameters, and this 
strongly suggests that the new gravity model must be a 
correct theory of gravitation. 

In addition, this model is applied to the description of 
the advance shift of the earth rotation around sun 0]. 
The predicted time shift for the earth rotation around 
the sun for one year is calculated to be 

(AT)t,, ~ 0.621 s/year 

which should be compared to the observed time shift in 
terms of the leap second 



(AT) 



exp 



0.63 ±0.02 s/year. 



The agreement between the prediction and the observa- 
tion is perfect, and this is discussed in detail in [J. 

It should be fair to make a comment on the general 
relativity [5]. However, the intrinsic problems of the gen- 
eral relativity are well described and discussed in detail 
in [ij, 1^. Therefore, we discuss only the numerical results 
of the wrong direction of the Mercury perihelion shift by 
the general relativity. The basic mistake arises from the 
calculation of the angular shift instead of the angular ve- 
locity shift. This is quite clear since the angular shift 
is not a physical observable. To understand intuitively 
a simple reason for the advance or retreat shifts of the 
Mercury perihelion, we discuss the shift in terms of the 
physical observable luT with to angular velocity and T 
the period. This quantity wT is basically proportional to 
the elliptic area of the Mercury orbit, and the new ad- 
ditional potential is repulsive and therefore it gives a bit 
larger area than the Newton solution, while the general 
relativity is attractive, and therefore it gives a smaller 
area. This indeed corresponds to the advance or retreat 
shifts of the Mercury perihelion, and we will show it later 
more in detail. 



II. NEW GRAVITY MODEL 
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In this section, we briefly review the new gravity model 
[HQ. First, we write the Lagrangian density 

£ = i-ipj^^d^ip — e-ip^^A^ij] — m(l + gG)^^^^ 



2 



(2.1) 



where and denote the electromagnetic field and 
gravitational field, respectively. The gravity is a mass- 
less scalar field which couples to the fermion field -0 in the 
mass term. This Lagrangian density is gauge invariant, 
and therefore it keeps the most important local symme- 
try. 



A. Photon-Gravity Interaction 

As shown in [H, Q , photon interacts with gravitational 
field at the fourth order of Feynman diagram. The de- 
tailed expression is given in [l|,Q, and therefore we write 
only the result of the equation of motion for photon 
A(t, r) under the gravitational field 



9t2 



where M, mt and a denote the mass of the gravitational 
center, the total mass of different kind of fermions and 
the fine structure constant, respectively. An interesting 
question may be as to whether there is some experiment 
which can measure the photon gravity interaction. We 
believe that, if one makes use of the satellites, there may 
well be some chance to measure it because of the laser 
technique. From one satellite to the other, one can shoot 
some laser beam with very low energy. In this case, this 
laser beam may well be scattered by the earth gravity, 
and the scattered laser may be observed at the other 
satellite detector. 



B. Proper Dirac Equation in Gravitational 
Potential 

For a long time, the understanding of the Dirac equa- 
tion in the gravitational potential has been incomplete. 
That is, a proper Dirac equation for fermions was not 
known in the presence of the gravitational potential. 
This was not a healthy situation for theoretical physics 
since, it indicates that, even the Newton equation with 
the gravitational field cannot be derived properly. It is, 
in this respect, quite important that we have now the 
Dirac equation for fermions in the gravitational poten- 
tial. From the Lagrange equation of eq.(2.1), we obtain 



iYdf.i' - ejf'A^tl) - m(l -I- gg)ij = 0. 



Also, we can write the equation of motion for the gravi- 
tational field 



(2.36) 



Under the static approximation for the gravitational 
field, we can derive the Dirac equation of eq.(l.l) for 
fermions in the gravitational potential. 



III. NON-RELATIVISTIC GRAVITATIONAL 
POTENTIAL 

A. Foldy-Wouthuysen transformation 

Now, the Hamiltonian of the Dirac equation in the 
gravitational field can be written as 



H — — zV ■ a + m — 



'''^^'-] p. (3.1) 



r J 



This Hamiltonian can be easily reduced to the non- 
relativistic equation of motion by making use of the 
Foldy-Wouthuysen transformation Q. Here, we only 
write the result in terms of the Hamiltonian H 

p2 GmM 1 GmM . 1 GMm , 



2m 



2irfl r 



(3.2) 

where the last term denotes the spin-orbit force, but we 
do not consider it here. Now, we want to make the clas- 
sical limit to derive the Newton equation. In this case, 
it is safe to assume the factorization ansatz for the third 
term, that is. 



1 GmM 
2m2 r 



-P = 



1 GmM 
r 



{p') ■ (3.3) 



By making use of the Virial theorem for the gravitational 
potential 



GmM 



(3.4) 



we obtain the new gravitational potential for the Newton 
equation 

GmM 1 (GmMV 

V{r) = + — ^ (3.5) 

r Zmc \ r / 

where we explicitly write the light velocity c in the last 
term of the equation. 



B. Large Component Treatment 

The new gravitational potential can also be obtained 
from the equation for the large component in the Dirac 



(2.3a) equation. By denoting the wave function ip as ip = 



we can rewrite the Dirac equation 
GmM 



X 



E 



cr-px = 0, 



, , GmM\ 
cr ■ pep - [ E + m ;^ = 0. 



(3.6) 
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By denoting E = m + Enr, we obtain the equation for 



p2 GmM 1 /GmM\ 1 GMm 



2m 



2m 



(EnJ)'^ 

2m 



— En 



(3.7) 



where, in the last term, we replace Enr by Enr which is 
the eigenvalue of the unperturbed Hamiltonian. In this 
case, we obtain the total gravitational potential 



V{r) = - 



GmM 1 /GmMV 
r 2mc^ \ ^ J 



which is just the same as eq.(3.5). As mentioned above. 
Coulomb potential has no such effect from the non- 
relativistic reduction, and therefore, this type of the ad- 
ditional potential is special for the scalar potential. 



IV. ANGULAR VELOCITY SHIFT OF 
MERCURY AND GPS SATELLITE 

The Newton equation with the new gravitational po- 
tential can be written as 

GmM P G^M^m , 

mr = — + - + (4.1) 

mr^ c r 

Therefore, we can introduce a new angular momentum L 
as 



(4.2) 



Further, we define the angular velocity oj and radius R 
by 



R 



GMm2(l -£2)i 



(4.3a) 



where e denotes the eccentricity. Correspondingly, we 
can define a new angular velocity U, associated with lj as 



G^M^ 



u;' 1 + 



G^M^ 

where 77 is defined as 

G2M2 



^'(1+^) 



C^R^LU^ 



(4.35) 



(4.4) 



The equation (4.1) can be immediately solved, and one 
finds the solution of the orbit 



A 



1 + ecos {■jip) 



(4.5) 



where A and e are given as 

£2 21JF 



A = 



GMm^ 



e=^ I 



m{GmMf 



(4.6) 



Physical observables can be obtained by integrating tp 
— T- over the period T 



It; r2-K 



— dt^ r^dip = A^ 

rnJo Jo Jq (1 + ecos (^93))^ 



1 



■dip. 



(4.7) 



This can be easily calculated to be 

LoT ^ 27r(l + 2r/) (1 - erj) ~ 27r{l + (2 - e)ry} (4.8) 

where e is assumed to be small. Therefore, the new grav- 
ity potential gives rise to the advance shift of the angular 
velocity, and it can be written as 



(2-£)77. 



(4.9a) 



th 



This is a physical observable which indeed can be com- 
pared to experiment. Here, it may be interesting to com- 
ment on eq.(4.8). It shows that one may also discuss the 
physics of the deviation in terms of the time shift. 



AT 



(2-£)77. 



(4.96) 



th 



In fact, the shift of the GPS satellite is often described 
in terms of the time shift. 



A. Mercury Perihelion Shift 

The Mercury perihelion advance shift A9 is well known 
to be 



AO ~ 42 " per 100 year. 



(4.10) 



Since Mercury has the 0.24 year period, it can amount 
to the shift ratio SO 

SOobs = (—] ^ 7.8 X 10-«. (4.11) 

The present theoretical calculation shows 
G2m2 



V 



,2 ^4,, ,2 



2.65 X 10" 



(4.12) 



where the following values are used for the Mercury case 
i? = 5.73x10^%!, A/ = 1.989x10^" kg, w = 8.30x10"^ 
Therefore, the theoretical shift ratio 69th becomes 
'Alu' 



4.8 X 10" 



(4.13) 



th 



which should be compared to the value in eq.(4.11). As 
can be seen, this agreement is indeed remarkable since 
there is no free parameter in our theoretical calculation. 
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GPS Satellite tj-Shift 



where 7 is found to be 



Many GPS satellites which are orbiting around the 
earth should be influenced rather heavily by the new 
gravitational potential. The GPS satellite shift can 
be estimated just in the same way as above, and we ob- 
tain 



1.69 X 10" 



(4.14) 



where we employ the following values for the GPS satel- 
lite i,ai 

R = 2.6561 X lO'^ m, M = 5.974 x lO'^^ kg, 



uj = 1.4544 X 10""* (4.15) 

since the satellite circulates twice per day. Therefore, the 
advance shift of the GPS satellite becomes 



Aw 

UJ 



3.4 X 10" 



(4.16) 



th 



This should be compared to the observed value of 



4.5 X 10" 



(4.17) 



exp 



As seen from the comparison between our calculation and 
the observed value, the present theory can indeed achieve 
a remarkable agreement with experiment. 



V. PREDICTION FROM GENERAL 
RELATIVITY 

Here, we briefly discuss the calculated result by the 
general relativity. For the Mercury perihelion shift, the 
result is quite well known, and it can be written in terms 
of the angular shift. In fact, the angular variable is 
modified by the general relativity to 



cos (f 



cos(l — 7)(^ 



(5.1) 



(5.2) 



This change of the shift in the angular variable explained 
the observed Mercury perihelion shift. However, as can 
be seen from eq.(4.5), this effect vanishes to zero in the 
case of £ = 0, that is, for the circular orbit. This is, of 
course, unphysical in that the effect of the general rela- 
tivity is valid only for the elliptic orbit case. In Newton 
dynamics, the angular momentum i is the only quantity 
which can be affected from the external effects like the 
general relativity or the additional potential. As can be 
seen from eq.(4.5), the effect of the angular variable in 
eq.(5.1) induces the change in A of eq.(4.5) as well, and 
this is indeed inevitable. Therefore, we should calculate 
the physical observables as to how the general relativity 
can induce the perihelion shift, and we find 



wT ~ 27r{l - 2(2 - e)7}. 



(5.3) 



Unfortunately, this is not the advance shift of the Mer- 
cury perihelion, and instead, it is a retreat shift. 

This is a reflection of the important consequence in 
physics that one should evaluate physical observables in 
order to compare theory to experiment. In this respect, 
the angular shift calculated by the general relativity [1, 0] 
is not a physical observable. In fact, as can be easily 
seen, the angular shift at one fixed point of the perihelion 
cannot take into account the effects that should arise 
from all the orbital trajectories of the Mercury since it 
is rotating around the sun. The shift must be the result 
of the integrated quantity over the orbit, and it is not 
a quantity just at one point of the trajectory. Since the 
observed shift of the Mercury perihelion was obtained 
in comparison with the one that was made some years 
before, there is no chance that the observed perihelion 
shift is directly related to the angular shift in eq.(5.1). 
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